r 

AO-A'UO  479  WISCONSIN  UNXV-MAplSON  MATHEMATICS  RESEARCH  CENTER 
ON  AN  AXISYMMETRIC  FREE  BOUNOARY  PROBLEM. (U) 

AUG  61  S  ZHOU 
UNCLASSIFIED  MRC-TSR-2270 


P/G  20/11 


DAAB29-80-C-0041 

NL 


"W03  3111  m 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 

610  Walnut  Street 
Madison,  Wisconsin  53706 

August  1981 


Received  April  8,  1981) 


Approval  for  public  release 
Distribution  unlimited 


Sponsored  by 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 


and  National  Science  Foundation 

Washington,  DC  20550 


82  02  03  059 


UNIVERSITY  OF  WISCONSIN  -  MADISON 
MATHEMATICS  RESEARCH  CENTER 

ON  AN  AXISYMMETRIC  FREE  BOUNDARY  PROBLEM 

Shu-Zi  Zhou* 

Technical  Summary  Report  #2270 
August  1 981 
ABSTRACT 


"The  axisymmetric  elastic-plastic  torsion  of  a  shaft  of  general  shape 
subject  to  the  Hencky  consistency  condition  with  the  von  Mises  yield  function 
is  considered.  It  is  proved  that  the  Haar-Kannan  principle  is  valid  in  this 
case,  and  that  the  problem  is  essentially  two-dimensional.  The  problem  is 
reformulated  as  a  variational  inequality,  and  the  existence  and  uniqueness  of 
the  solution  is  studied. 
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SIGNIFICANCE  AND  EXPLANATION 


When  a  shaft  of  circular  cross  section  is  subjected  to  a  torque,  a 
istic  enclave  may  appear.  The  boundary  between  the  elastic  region  and  the 
plastic  region  is  unknown.  It  is  a  so-called  free  boundry  problem. 

Assuming  that  Hencky's  consistency  condition  with  the  von  Mises  yield 
function  is  satisfied,  we  can  prove  that  the  Haar-Karman  principle  is  valid, 
which  means  that  the  strain  energy  must  be  minimized  subject  to  the  constraint 
that  the  stress  should  not  exceed  its  permissible  limit:  •  We  show  that  the 
problem  is  essentially  two-dimensional,  and  give  two  kinds  of  variational 
formulations  of  the  problem:  one  for  the  stress  field,  the  other  for  the 
stress  function.  The  existence  and  uniqueness  of  the  solution  of  the 
variational  inequality  for  the  stress  function  is  proved. 


The  rei" jnsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sumr.ary  lies  with  M*C,  and  not  with  the  author  of  this  report. 


ON  AN  AXISYMMETRIC  FREE  BOUNDARY  PROBLEM 


Shu-Zi  Zhou* 

1 .  Introduction 

The  elastic-plastic  torsion  of  shafts  is  one  of  the  classic  free  boundary 
problems.  For  the  case  of  that  the  cross-section  is  constant,  it  has  been 
studied  deeply  by  using  variational  inequalities  during  the  last  decade  'see , 
for  instance.  Ting  [1971,  1976],  Bresis  and  Sibony  [1971],  Lanchon  [1974], 
Friednan  [I960],  Pozzi  [1980]).  Recently,  Cryer  [I960]  has  considered  the 
case  in  which  the  shaft  has  variable  cross-section  and  rotational  syssietry. 

He  has  proved  the  existence,  uniqueness  and  regularity  of  the  solution  of  the 
variational  inequality  problem  for  the  stress  function  under  some  assumptions . 
He  has  assumed  that  the  function  it  describes  the  generator  of  the  rotational 
shaft  is  monotone.  In  this  paper  we  consider  the  case  in  which  the  generator 
may  have  more  general  shape;  give  two  kinds  of  variational  formulation  of  the 
problem:  one  for  the  stress  field,  the  other  for  the  stress  function;  prove 
that  Haar-Karman  principle  is  valid  and  the  problem  is  essentially  two- 
dimensional  under  the  so-called  Hencky's  conditions;  study  the  existence  and 
uniqueness  of  the  solution. 

I  am  grateful  to  Professor  C.  Cryer  for  many  valuable  discussions. 
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2 .  Classical  problem 


Let  Figure  1  represent  a  shaft  of  circular  cross  section  with  equal  and 
opposite  pure  torques  T  applied  at  the  ends .  Assume  that  the  material  is 
homogeneous,  isotropic,  and  elastic-perfectly  plastic,  that  there  are  no  body 
forces ,  and  that  there  are  no  external  tractions  on  the  lateral  surface .  Our 
aim  is  to  find  the  resulting  stress  distribution. 


Figure  1 

Setting  up  a  cylindrical  coordinate  system  in  which  the  z-axis  coincides 
with  the  center  line  of  the  shaft  and  the  origin  lies  in  a  end  of  the  shaft, 
we  assume  that  the  equation  for  the  generator  is  r  R(z)  .  Let  (see  Figure 
2) 

ft  =*  {(z,r):  0<z<L,  0<r<  R(z)} 
rQ  =  { (z ,r ) :  0  <  z  <  L,  r  -  0} 

T1  «.  { (z  ,r) :  0  <  z  <  L,  r  -  R(z)} 

T21  *»  {(z,r)s  0  <  r  <  R(0 ) ,  z  -  O) 

”  f(z,r):  0  <  r  <  R(L) ,  z  -  l) 


-2  ■ 


7^* 


t 


Figure  2 


Then  we  have  (Eddy  and  Shaw  (1949]). 

1  —  2  „ 

Classical  Problem  i  Find  stress  function  v  e  c  (0)  n  c  (w)  such  that 

2 

|Vv|  <  kr  in  ft  (k  is  a  constant  given) 

.  _  3  -3  3v.  3  -3  9v, 

Av  =  "  Si  (r  3l>  "  3?  (r  3?’  “  0 

in  n  ^  {(r,s)t  |?v|  <  kr2} 
v  *  0  on  T  ,  v  *  T/2W  on 

3v  n  P 

^  -  °  on  r2  . 

This  formulation  of  the  problem  is  based  on  the  von  Mises  yield  criteria, 
and  on  the  so-called  semi-inverse  method  given  by  Saint-Venant  (see,  for 
instance,  Timoshenko  et  al.  (1951,  p.  259,  p.  306})  in  which  one  assumes  a 
priori  that  there  are  no  radial  and  axial  displacement.  This  fact  will  be 
proved  in  section  4  of  this  paper. 

j .  Haar-Karman  principle 

The  argument  we  use  in  these  two  sections  is  similar  to  that  in  Lanchon 
(1974]  for  the  case  of  constant  cross  section. 

From  now  on  we  assume  that 

R(z)  is  piecewisely,  continuously  differentiable 

R’(0>  >*  R'(L)  t  V  (3.1) 


which  implies  that  ft  is  strongly  Lipschitz  domain.  Then  the  three 
dimensional  region  ft  occupied  by  the  shaft  is  also  strongly  Lipschitz 
domain .  Its  boundary  is 


an*  -  r;  u  r*t u  r;2 

*  *  ♦ 
where  T  is  the  lateral  surface  while  and  " ^  are  the  end  surfaces  • 

Denote  by  ur#  uQ  and  uz  the  components  of  the  displacements  in  the 

radial,  tangential  and  axial  directions  respectively*  Let 

u  -  tur,u0,uziT,  0  -  fVVVaez'ar*'arelT  ' 

e  -  [e  ,efl,e  ,e  ,e  ,e  fl]T 
r  8'  z  8z  rz  r0 

where  0-stress  field,  e-strain  field.  Then  we  have  (see,  for  instance, 
Timoshenko  et  al.  [1951,  pp.  305-308]) 


3u  u  3u0  3u 

Er  “  TT'  e0  “  71  +  r  W'  EzmTT 


3u0  1  3u 


3u  3u 
r  .  z 


3u  3u.  u. 

1  r  8  8 


e0z  *  JT  +  7  50“'  erz  "  7T  +  5T'  ere  “  7  W  *  SF"  '  F* 


9a 
_ r 

TT 


.  3a 
1  r8 

7TT 
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in  ft 


(3.2) 
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3O0 

36 
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,  i  3aez 
*  7  ~7T 


o 


in  ft 


(3.3) 


vie  have  the  boundary  condition  as  follows i 

an  +o  n  "0 
r  r  rz  z 

a  0n  +  aft  n  -  0 
r8  r  0z  z 

_* 

a  n  +  a  n  *  0  on*,  (no  external  tractions 
rz  r  z  z  1 

°rz  *  °0z  “  °z  =  0  °n  r21  U  r22  <pure  torque) 

/  .  rO0zda  “  f  •  ^Oz**  '  T 
21  22 


) 


(3.4) 

(3.5) 
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(3*6) 


where  (nr,  ng,  nK)  -  n  ia  the  outer  normal  of  ^0  .  (obviously,  n  is 

*  * 

well-defined  a.e.  on  3(J  because  of  (3.1),  and  n^  -  0  on  ««  ,  nr  *  0  on 

T*  U  r*  ) . 

21  22 

Assume  that  the  resulting  stress  field  o°  satisfies  the  Hencky's 
consistency  conditions  (Lanchon  (19743) 

P(o°)  <  0 

0  0  , 
e  -  AO  +  X 


lT(o-o°)  <0,  V  a  e  M1 


(3.7) 

i 

where  F(o)  is  the  von  Mises  yield  function 

•>  •> 

(3.8) 

X  a  [X  , •  •  • , X  lr,  m  *  {o  s  F(O)  <  o),  and  A  is  the  matrix  in  the  Hooke's 
r  rs  1 

law  (Timoshenko  et  al.  [1951,  p.  7,  p.  66}) 


F(o)  -  \  (o2  ♦  o*z  *  a2cz  +  a2fl)  -  1  (or  +  o9  +  oz)2  -  k2 


1  -V 
-V  1 
-V  -V 


-V 

-V 

1 


2(1+V) 

(2(1 -tv) 

2(1+v). 


(3.9) 


where  E  is  the  modulus  of  elasticity  and  v  is  the  Poisson's  ratio  with 

0  <  v  <  1/2. 

Let 

M2  -  {o  e  M1  s  o  e  tH 1 (0*) ] 6,  0  satiofies  (3.3)  -  (3.6)} 
where  the  o  in  (3.4)  -  (3.6)  is  the  trace  of  a  on  98  (see,  for 

instance,  Necas  [1967,  p.  15]). 

Proposition  3.1.  If  (3.2)  and  (3.7)  are  valid  for  o°  e  Mj,  then  the 
Haar-Karman  principle  is  valid  in  this  case,  that  is 


5- 


where 


J  (0  )  -  min  J  (o) 

o&*2 


J.  10)  -  ^  *  °TA<*dV  . 

I  4.  n 


(3.10) 


Proof :  v  ®  6  Mji  we  have  (noting  that  A  is  symmetric) 


0  ^  t  T  U  T  0 

-  J.(<J  )  -  j  /  *(«  AO  -  (O  )  AO  )  <JV 

n 

1/  0T  0  »  0  I  0 

-  4  J  ,<0-cVa(<J-0U)<JV  +  /  #<0VA(0-0U)dV  . 

2  n  n 

It  follows  from  Gerschgorin's  theorem  (Varga  [1962,  p.  16])  that  all  of  the 
eigenvalues  of  A  are  not  less  than  (1-2V)/e  >0.  Hence 

xTAx  >  (1-2V)x2/E,  V  x  e  R6  (3.11) 

and  we  have 


J AO)  -  J.(o°)  >  /  #  (o°)TA(o-o°)dV 
1  1  fl 


r  ,  (E  )T 

(a-o° 

)dV 

n 

3u° 

0 

u 

(ar  - 

o°)  +  ( 

—  +  - 

*r  j 

»u° 

3u° 

3u°. 

33^  (°9x 

"°9z) 

+  (3 T 

+  3T' 

(since  (3.7)) 


3u 


.  0  .  0  0 

«U  OUg  Ug 

+  - — J  (Org-c“g)]rdrd9dz  (since  (3.2)) 


/  «,(u°((o  -o°)n  +  (o  -o°  )n  )  +  uS{(®  fl-a°)n  +  (aA  -a°  )n  ) 

gjj*  r  r  r  r  rz  rz  z  9  r9  r9  r  0z  9z  z 


0  0  0  ,  0  3  0 

+  uz((°rz~°rz)nr  +  (VaZ)nz))dS  _  (ur(37  (W  + 

1  3  0  3  0  0.3  .  0  . 

7  35  (Or0‘°r9)  +  37  (0rz_°rz>  +  - r - ^  +  u0(37  (°r0-°r9)  + 


7  Is  <v!>  *  h  *  .M?  <«„<)  * 


7  Is  "’ei"°8,)  *  51  lV“!1  *  7  ror”“1*> 


«  0  .  (since  (3.3)  -  (3.5)) 

Q  «E  «D  • 

About  more  general  discussion  on  the  Haar-Karman  principle,  see  Martin 
(1975,  pp.  733-736]. 


4 .  Variational  formulation  of  the  problem 


Proposition  3.1  suggests  the  following  variational  formulation  of  our 


problem . 


Problem  (A)  Find  o  e  M2  such  that 


J  («)  -  min  J^d) 

oen2 


where  J1  ie  defined  by  (3.10). 

It  is  easy  to  show  the  following  assertion. 

Proposition  4.1.  Problem  (A)  is  equivalent  to  the  variational  inequality 

o°  e 


a,**0,  a-°°)  >  0 


v  a  e 


where 


»,(o1,o2)  -  /  (o1 )TA  o2  dV 

1  n 


Proposition  4.2.  Problem  (A)  has  at  most  one  solution. 


Proof:  If  0°  and  o1  are  solutions,  then  we  have  by  (4.1) 


a^o0,  o1-*0)  >  0 
a.  io1 ,  o°-a1 )  >  0 


-T 


By  virtue  of  (3.11)  we  obtain 


,<0e.'*>Q.rl  "  'h  ^  *  ff8z*  dV| 


ft 

<  JL  Ia?  I  _  *1*1  „  ^  -  -_L  la?  I  ^  •I'fi  „ 

21T  6Z  L2(fl  )  L2(fl  )  737  0Z  L2(fl  )  r  L2(fl) 

* 

Hence  0^  (more  precisely ,  its  restriction  in  ft)  belongs  to  the  dual  space 


of  IT  (ft)  and  then 
r 


°8z  e  Lr(fl)  ‘ 


So 


/  *  (o0z)2<iv  ■  C  d6  /fl  <o0z)2rdrdz  -  2z,0Sz,22  „  < 


L*{0) 


•  2  * 

i .e .  e  L  (ft  ) .  Similarly,  we  have 


3o“ 


**3iC- ir*  •Ml  ■ 


0  * 

<  *  k 


Given  v?  e  C0<ft),  we  have 


if  3oe 

<v'  *>a,r  “  k  C  s r  d0)*rdrdz 


/  -J—  f  dV  *  •==■  /?"  d9  /_  — y—  ^rdrdz 

2  w  o*  ®r  2*  0  ft  ®r 


1  r21f  f  0  w»r)  .  . 

—  J  .  do  Jn  o  — - -  drdz 

2»  0  ^  ft  “z  or 


0  3  (<Pr) 


1  r  ,r2»  J>  ,a.  3<^r-)  .  . 

-  2*  ]n  <Jo  °0z  d6)  T?  “ 

* 

-  k  °\z  ^r1  “  “  ^  it  ”rdrdz 


3qez 


It  means  that 


<39  e  L2  (ft* ) 
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Similarly,  we  obtain 


3  at  3o? 

Sz  1  r2*  9z  ..  a  2<0* 

T7“  *  2*  fo  “57  <je  e  L  (£1  )  . 


But 


•V.  , 

"55"  *  0  • 

*  1  *  1  * 

Hence  Og^  e  H  (»  ) .  The  same  argument  indicates  that  o^g  e  H  (w  )  .  and  we 


obtain 


*  1  *  e 

o  e  IH  (d  )] 


Turn  to  (3.3)  -  (3.6).  Check  the  second  equations  of  (3.3)  and  (3.6)  for 


* 

O  : 


3o*  1  3ofl  3afi  20*  3a°  3o°  20° 

_J*  +  1  »  +  °£  +  _£°  .  ’  /f  (  *£  +  °£  +  — Ei)  d6 

3r  r  do  3z  r  2*  '0  v  3r  3z  r 


3o" 


1  r2»  1  0 

2?J0  7srd0 


(since  (3.3)) 


(periodicity  of  a  in  6) 


V 

22 


/ .  r<^  /r  -t d9)  rd0,dr 

22 


57  /f  /;,L>  *  r2  d6 

2»  2  0  2  0  0  0Z 


(R(L)  f2*  2  0  .  r  0  _ 

‘o  *  ‘o  *  °0Z  d0  *  ro0z  * 

22 


The  rest  of  (3.3)  -  (3.6)  is  clear. 

Now  we  prove 

•  *7  *  5  5 

f(o  )  -  (°er)  +  «>re>  -  *  <  o 


(4.4) 


a>  -  a 


0.  6. 


+  ®  ft 


Me  have 


1  f2ir  ,  0  *  0  * 

=  —  I  a.  a a  +  a  „a  a)d9 
2ir  '0  8z  9z  r9  r0 

<  h  (C<(09z)2  +  (Ore)2)d0)1/2  *  ^  <(°8z>2  +  <°°0)2><i6)1/2 

“  7^  t(°9z)2  +  (0r9)25,/2‘  </T  t(«°9z)2  *  (oJe>2lde)V2 

(O0*z)2  +  ,Or0)2  <  h  ft  f(°2z)2  +  ^6>2Id0  * 

On  the  other  hand,  we  have 


i  «e  • 


(°9z)2  +  (0r9)2  ^  F(°0)  +  **  *  *2  • 


Therefore,  we  obtain 


«v2  ♦  ^;9>2  <  h  % - *2 


i.e.,  <4. 4)  is  valid,  and  a  e  M^. 

0  * 

Finally,  we  compute  J^o  )  -  J ^(0  ).  Clearly, 

(oVao0  -  ^  [(O0)2  +  <0°)2  +  (0°)2  -  2V(0°0°  +  oV  +  o°c°)] 

+  <)2  ♦  «£/»  >  ^  <<>2  ♦  «>v-> 


It  follows  from  (4.5)  that 

*2  .  *2 


(since  0  <  v  <  — )  , 


/  .  I +  ^e*2^  «  /  *  t <°L>2  ♦  • 


9z 


r9 


9z 


r9' 


.K.nce 


0  *  r  0T0  *  T  * 

2(J,{0  )  -  J  (<J  ))  -  J  #  [<oVaou  -  (0  )  AO  JdV 

I  S  pk 


2  ( 1  tv ) 
E 


</  •  [(°9z)2  +  (Or0)2]dV  '  ^*t(°9z)2  +  Kg)2^  >  0 


a  n 

(4.3)  has  been  proved,  and  the  proof  of  the  theorem  is  completed. 


•11- 


(4.5) 


Q.E.D 


The  basic  idea  of  the  above  proof  is  the  same  as  that  of  the  so-called 


semi-inverse  method. 

This  theorem  enables  us  to  take  the  set 

H  *  {d  8  H,  :  3  =a„-0  =  0 

2  r  6  z  rz 


0  in  Q  } 


as  the  set  of  the  adnissible  stress  vectors  instead  of  M2  in  problem  (A)  . 

Then  we  have  the  following  problem. 

Problem  (B) .  Find  e  N  such  that 

J  (o°)  -  min  J^o)  . 

aeS 

Obviously,  problem  (B)  has  at  most  one  solution;  if  o®  is  the  solution 
of  problem  (A),  then  it  is  also  the  solution  of  problem  (B)  . 

Remark  4.1  .  If  a  e  S,  then  it  follow  from  (3.3)  that 


3o 


Therefore, 


where 


'8z  r0 

T§~  *  ~7B~ 


V®>  ■  4  *  "s,11  drdz 


5. 


The  variational  problem  for  stress  function 
If  problem  (B)  has  solution  e  [C°(fl  ) ] 6 ,  then  we  have  by  remark  4.1 

3°  e  n  ,  j  (a°>  “  Bin  J0(a>  (c> 

®eN1 


where  N.,  is  a  subset  of  N: 

n  -  {o  e  n  ;  0  e  [c°(n*))6} 

V  o  e  N  ,  by  virtue  of  (3.3)  we  have 

3  2  3  2 

7?  (r  °re)  +  71  (r  °9z)  *  0  * 

Then  it  is  easy  to  show  that  there  exists  v  e  H*(ft)  c'(fl)  such  that 
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in  n 


3v  2  3v  2_ 

W  "  r  °0z'  75 - -  °r9 


Hence 


•jj—  ■  0  on  Tp,  and  v  ■  on  Let  v  »  v*  -  C^.  Then 


3v  2  3v 

17  “  r  °9z '  -57 

V  ■  0 


■  -r  o  in  U 


rO 


on 


(5.1) 

(5.2) 


J0(O)  *  4  r~3[(!i:)2  +  (•|j)2]<Sr<i*  5  J(v>  • 

It  follows  from  (3.4)  and  (5.1)  that 

dv  3v  ,  .  ,  3v  ,  .  _ 

—  -  -  oos  (n,r)  +  -jj-  cos  (n,«)  -  0  on 

where  s  is  the  tangential  direction  of  1*^.  Therefore, 

v  -  C_  on  r 


v  ■  T/2w  on  T 


By  (3.5)  and  (5.1)  we  have 

0 

on 

It  follow  from  F (<X)  <  0  that 


2  * 


10  1 2  4,  2  2  .  c  J2.  4 

(Vv(  *  r  (dr0  +  aQz)  <  k  r 


|Vv|  <  kr2 


in  8 


Since  o  e  [H1^  )]®,  we  have 


0r9'  °0z  6  Hr(0) 


where 


H^(fl)  »  {v  0  L2(8)  s  e  L^(O)} 


(5.3) 


(5.4) 


2  “"1  * 

Since  (3.6),  (S.1)  and  (5.2),  we  have 

T  -  /**  d9  /*(0)  r2a0j!drd0  -  2w  /*(0)  dr  -  2*  v(0,R(0))  ,  (5.5) 

Combine  (5.4)  and  (5.5)  we  obtain 


(5.6) 


(5.7) 


(5.8) 


Hence 


L*  (ft)  is  defined  by  (4.3). 


-2  3v  -2  3v  .  „1,n. 
r  TTZ '  r  TT  8  H-<n>  » 


z  r 


(5.9) 


i 
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particular  ly 


(5.10) 


fa  r"3  /q  *_3<jf>2*<*«  <  “ 


~3/3vj2. 
z 


Finally  by  (5.1)  and  (5.2)  we  have 


V2  -  til  rV<*>2  <  /'  r4dr  .  /!  a2_dr  < 


Cr~ 


so 


Iq  r  5v2drdz  <  *  . 


(5.11) 


Now  let 


N2  -  (v  :  V  e  H2(fl)  nc1)!!),  (5.2)  and  (5.6)  -  (5.9)  are  valid}  . 

Then  it  is  easy  to  see  that  (5.1)  -  (5.2)  defines  a  biunivocal  map  from  N1 
onto  N2 ,  and  problem  (C)  is  equivalent  to  the  following  problem: 

(D) 


where  (noting  (5.3)) 


Vq  e  N2,  j<vq)  *  min  J(v) 
veN2 


J(v)  “  /Q  Pt(|^)2  +  (|j)2]drds 
.-3 


(5.12) 


P  -  r  -  .  (5.13) 

Now  we  enlarge  the  set  of  admissible  functions  of  the  variational  problem 
(D)  for  solving  the  problem  on  the  existence  of  the  solution.  Noting  (5.11) 
and  that  only  the  derivatives  of  first  order  appear  in  the  functional  J(v), 
we  Introduce  a  set  as  follows 

K  -  {v  :  v  e  H’(fl),  v  -  T/2W  on  1^,  |Vv|  <  kr2  in  «}  (5.14) 


where 


Hp(fl)  -  (v  e  L2(JJ)  :  e  L2(«)} 

L  (0)  »  {v  :  v  measurable,  Ivl  .  <  +“} 

Lp<n) 


(5.15) 

(5.16) 


with  norm,  respectively. 


Ivl  -  (/0  Pv2drdz]^2 

■>> 


Ivl  .  (Ivl2  *  ,3v-2 

Hp(fl)  Lp(0)  vt  Lp(SJ) 


,2-n.  oz  _  2 


.  ii:i2,  >h  , 

l2(H) 


-14. 


and  v  -  T/2*  is  in  the  sense  of  trace  (see,  for  instance,  Necas  [1967,  p 


15]).  K  is  a  natural  extension  of  Nj  because  of  that  according  to  theorem 
2.2  in  Cryer  [1980],  if  u  6  H^(Q),  then  (5.2)  (in  the  sense  of  trace)  and 
(5.11)  are  valid.  Therefore,  we  have  the  variational  problem  for  stress 
function  as  follows  . 


Problem  (B)  .  Find  vQ  0  X  such  that 

J(v  )  ■  min  J(v) 
0  v8K 

and,  equivalently. 

Problem  (F)  .  Find  Vq  such  that 


fvo8K 

a ( vq  ,  v-vQ)  >0,  ¥  v  e  K 


where 

,  3v*  3v*  3v* Jv" 

a(V  ,v")  -  Ja  ♦  £  .j2_)drd* 

Similarly  to  the  case  of  constant  cross-section,  we  introduce  the 

obstacle  problems  relevant  to  problesi  (F) .  There  are  two  obstacle  problems  to 

be  considered. 

Problem  (FI).  Find  v^  such  that 

V,  e  x. 


1 


a(v1#  v-v^  >  0 


v  v  e  x 


where 

K1  *  (v  8  Hp[8>  «  v  ■  T/2»  on  1^,  v  >  ^  in  H)  , 
and  is  the  solution  of  the  Cauchy  problem 
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^  «c  w 


I 12  -  k2r4,  *  <  T/2ir  in  A 


(5.17) 


♦  -  T/2* 


on  . 


Problem  (F2).  Find  v2  such  that 


v„  e  K, 
2  2 


a(v2'  v"V  *  °  V  v  6  K2 


-  (v  e  Hp($l)  tv*  T/2*  on  ,  v  <  t2  in  ft)  , 


and  is  the  solution  of  the  Cauchy  problem 


1»2  e  C2(fl) 


|V*2|2  -  *2r\  i|»2  >  0  in  0 


(5.18) 


*2  “  0 


°n  ro  . 


The  problem  (FI)  is  just  the  problem  (4.7)  in  Cryer  [1980],  there  the  solution 
for  (5.17)  is  also  discussed. 


6.  -  solution  of  the  Cauchy  problem  (5.18) 


Assume  *  is  the  solution  of  the  Cauchy  problem  (5.18).  Let  p  -  y— 

a+2  2 

q  *  .  Then  the  equation  is 

_  _  2  .  2  2  4 

F  =  p  +  q  -  k  t  *0  . 

He  have  along  the  characteristics  parametrized  by  s  in  (Courant  and 

Hilbert  [1962,  p.  78]) 

i6-’> 


(6.2) 


dr  3f 

dl-JS"  * 


3F  .  3F 


2  4 


aJT  ”  P  ■Jp  q  ■  2(P  +q)“2kr 

iE  .  _(p  I JL.  |r,  .  o 

to  p  3*  +  3z' 


dq  ,  3F  a  3f.  -.2  3 

to  “  ~(q  3*^  +  3?’  “  *  r  * 


(6.3) 

(6.4) 

(6.5) 


Since  4>2  e  C2(fi),  (6.1)  -  (6.5)  are  alao  valid  on  r  •  Given  (z^.O)  *  rQ. 

consider  the  characteristic  passing  this  point.  Let  the  parameter  value  of 
this  point  be  s  •  0.  Mb  have  initial  conditions  (since  1>2  »  0  on  and 


F  -  0) 


z(0)  ■  zQ,  r(0 )  -  p(0 )  -  q(0 )  “  ^(0)  *  0 


It  follows  from  (6.4)  that  p(s)  =0.  Hence  a  5  a  (since  (6.1)).  and 
q  ”  kr2  (since  F  >  0  and  <P2  >  0).  Therefore,  we  obtain  by  integrating 
along  the  characteristic 

d$  3$  3^ 

*2 <a!0  'r*  ”  "  ft  'ZT  *  ■  ^  *  +  a~  ^Ito 


'0  to 


0  l3a  to  3r  to 


*  /*  q  to  -  /g  q  dr  -  /J  kr2dr  -  kr3/3  . 

Since  (zq.0)  e  rQ  is  arbitrary,  we  obtain  the  solution 

♦2(z,r)  »  kr3/3  in  T5  . 


(6.6) 


7 .  Properties  of  the  set  K 

0,1  ys 

Denote  by  C  (»)  the  set  of  functions  which  satisfy  Lipschitz 
conditions,  that  is,  if  v  e  c0,1(fi),  then 


|v(p,)  -  v(p2) | 


h .(v.S)  5  sup  - 4 - - 

1  Pl.P2en  pi  p; 


(7.1) 


PA2 
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where  p,  -  (z^r,),  p2  -  (z2,r2),  and  |p1  -  p2l  - 
[(z,  -  z2)2  +  (r1  -  r2)2]^  . 

Lemma  7.1.  The  following  statements  are  equivalent: 

(a)  v  e  C0'1  (R) 

1  A  3v  3 v  ®  _ 

(b)  v  e  H  (0),  and  yp,  ^  @  L  (ft)  • 

Proof:  It  has  been  shown  (cf.  Adams  [1975,  pp*  i«s )  that  if  v,  y^, 

3v  •  _ 

^eL  (ft)  then 

h  (v,fl)  <  c  (Ivl  „  +  ly-l  «,  *  )  .  (7.2) 

1  1  l  (Q)  3r  L  *»>  H  I.  (0) 

dv  3v  00  i 

By  similar  argument  we  can  prove  that  if  t— ,  -r-  e  L  (fl) ,  and  v  e  L  (ft) 

or  vz 

then 

h.(v,fl)  <  C  ( *"5“l  a  .  )  •  (7.3) 

1  2  3r  h  (ft)  3*  L  (ft) 


Therefore,  it  is  clear  that  (b)  Implies  (a). 

0,1  . 

Now  assume  v  e  C  (**) .  Then  v  possesses  a  total  differential  a.e.  in 
ft  (Morrey  [1966,  p.  65]),  i .e . 

<\  a 

v(Pl)  -  v(p2)  -  [y^Hz,  -  *2>  +  [•^Hr1  -  r2)  +0  (|p1  -  p2l)  (7.4) 

a  .e .  in  ft 

3  v  3  v 

where  [y^]  ,  [y^]  are  the  partial  derivatives  in  the  usual  sense.  Clearly, 
they  are  measurable  in  ft.  It  follows  from  (7.1)  and  (7.4)  that 

I [y“l  I ,  lly-j^l  <  h^v/H)  a.e.  in  ft  .  (7.5) 

Hence  v  e  H1(ft)  and  (Morrey  [1966,  p.  63]) 

,3v  3v  ,3v,  3v  .. 

l77)  “  “  J?  •  <7-6> 


3v  3v 


by  (7.5)  we  obtain  y^,  y^  e  L  (ft),  i .e .  (b)  is  valid. 


Q  «E  .D . 


Proposition  7.2.  K  C  Kj  O  K2 . 
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Proof :  Given  v  0  K,  it  is  enough  to  prove  that 

♦,  <  v  <  *2  •  (7.7) 

By  Lemma  7.1  we  have  v  0  C^'^$).  Then  (7.6)  is  valid.  Since  v  6  H^O), 
is  absolutely  continuous  in  r  on  0  <  r  <  R(s)  for  almost  all  value  z 
(Morrey  [1966,  p.  66)).  Therefore,  noting  that  v  ■  0  on  T  t  we  obtain 

v(r,z)  <  |v(r,z)|  »  l/*  t|^Jdr|  ■  l/g  |“  dr| 

<  /5  |Vv|dr  <  /Jj  kr2dr  “  **«•  ®  * 

The  seoond  part  of  (7.7)  has  been  proved.  Mow  prove  the  first  part.  The 
system  of  the  characteristic  equations  for  the  Cauchy  problem  (5.17)  has  the 
same  form  as  (6.1)  -  (6.5).  Let  the  parameter  value  B  -  0  correspond  to  the 
point  on  r  .  Then  it  follows  from  (6.2)  that  the  point  in  0  corresponds  to 
the  negative  value  of  s.  Noting  that  v  has  a  total  differential  a.e.  in 
(1  and  that  (7.6)  is  valid,  we  have  along  the  characteristic 

v(z(0)  ,r(0) )  -  v(z(s)  ,r(s) )  -  /J  ~  -  /J  (|f  2  +  If  2)d“ 

-  l°a  (|^  2p  +  |^  2q) ds  <  2  !°a  |Vvl*(p2  +  q2)^2  ds 

<  2  /°  k2r4ds  •  a*  ■  ^(zJOKrtO))  -  1^  (z(s)  ,r(s>  >  , 


But  v(z(0)  ,r(0) )  -  *  (z(0),r(0))  -  T/2*.  Hence  v  >  ^  . 


Q  eE  eD  e 

Proposition  7.3.  The  sets  K,  X1  and  K2  ar®  convex  subsets  in 

"J(°)* 

Proof ;  Given  v,  w  6  K  and  A  with  0  <  A  <  1,  we  have 


|V[Av  +  (1-A)w)  |  -  [  (A  +  (1-A)  |^)2  +  (A  |^  +  (1-A)  Ij)2]^2 
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<  ((X  | v-)2  +  a  |^)2i1/2  +  ui-a)2(-|2)2  +  n-x)2(|j)2]1/2 


-  X|Vv|  +  ( 1 -X) |?w|  <  kr  . 

By  the  linearity  of  the  trace  operator  (Necas  [1967,  p.  IS})  we  have 

Xv  +  (1-X)w  -  Xt/2*  +  ( 1-X)T/2*  -  T/2*  on  I*  . 

Hence  Xv  +  (1-X)w  e  K  and  K  is  convex. 

Let  (v^)  bo  a  Cauchy  sequence  in  X .  Since  8^(0)  is  a  Banach  space 

(Cryer  [19601  ),  there  exists  v  e  H  "*  ftt)  such  that  v  ♦  v  in  hVo).  It 

P  n  p 

follows  from  the  continuity  of  traoe  operator  (Necas  [1967,  p.  15])  that 
v  -  T/2*  on  r  .  By  well-known  subsequence  argument  we  obtain  that 
|Vv|  <  kr2  in  0.  Then  X  is  closed. 

The  aonclusion  about  and  Kj  can  be  proved  by  similar  argument. 

Q  .E  .D  . 


8 .  Solution  of  the  problem  (F) 

At  first  we  solve  the  problem  (F2 )  it  suggests  the  solution  of  the 

problem  (F)  .  Let  (cf.  Cryer  [1980,  Remark  5.3]) 

R  “  min  R(z)  (8. 

0<z<L 

kQ  -  (3T/2»)R  “3  .  (8. 

We  need  a  lemma  it  may  easily  be  shown  by  a  well-know  theorem  (Adams 
[1975,  p.  54]). 

Lemma  8.1 .  If  u  e  C°(0)  n  H^fO),  then 

tr  u  -  u  on  38 

where  tr  u  is  the  traoe  of  u  on  90. 

Proposition  8.2.  If  k  <  kg  then  the  problem  (F2)  has  no  solution.  If 
k  >  kQ  then  it  has  a  unique  solution. 


.  y ' 1  • 


Proof  i  There  exists  a  z  such  that 


0  <  r  <  L,  R(z)  -  R  . 

If  k  <  kg  then  by  (6.6)  and  (8.2)  we  have 

if>2(z,R)  -  kR3/3  <  K0i3/3  -  T/2ir  .  (8.3) 

Thus,  there  exists  a  real  number  c  and  an  open  neighborhood  £  of  the  point 


(z,R)  such  that 


♦2  <  c  <  T/2*  in  E1  5  Z  n  fl  . 


Assume  that  K2  is  nonempty.  Take  v  8  Kj.  Then 

v  <  c  <  T/2*  in  E  , 

and  there  exists  a  sequence  (vice  (0)  such  that 

n 


lv  -  vl 
**  !T(fl) 


0  (n  ♦  -)  . 


From  the  construction  of  vR  (Adams  (1 975,  pp.  54-56])  we  see  that  there 
exists  Ez  C  and  a  real  number  C  >  C  such  that  mean  T  >  o,  where 
T  m  n  T  ,  *nd  that 

v  <  C*  <  T/2*  in  T„  .  (8.5) 

n  2 

Then  by  the  continuity  of  trace  operator  and  (8.4)  we  have 

lv  -  tr  vl  <  lv  -  tr  vl  ,r  .  <  C.lv  -  vl  ,  ♦  0  . 

"  L2(T  )  n  W  1  "  H1 (Q) 

Hence  there  exists  a  subsequence  {v^}  which  converge  to  tr  v  a.e.  on  T  , 


and  by  (8.5)  we  obtain 


•  • 

tr  v  <  C  <  T/2*  on  T 


This  contradicts  that  tr  v  ■  T/2*  on  r  •  Therefore,  K2  empty,  and  the 
problem  (F2)  has  no  solution. 


If  k  >  kg,  then  let 


v  *  min  (ijij,  T/2*)  in  (5  . 


that  v  e  K2.  Clearly,  v  <  i|»2,  and 


v  -  kr  /3 


v  -  T/2 * 


in  £3  n  (r  <  d) 
in  fl  n  (r  >  d/ 


'  .  v 


I 


where  <3  -  (3T/2*  k) 


So  we  have  (Gilbarg  and  Trudinger  [1977,  p.  145]) 


1/3 


v  e  H^ft), 


3v 

3* 


0 


in  ft 


gp  -  kr^  in  ft  n  (r  <  d) 
gp  ■■  0  in  ft  n  (r  >  d) 


(8.7) 


It  is  easy  to  see  that 

«v«  <  *  . 

Hp(fl) 

Hence  v  e  Hp(ft).  Clearly,  v  e  C^(^)  ^  H^(ft),  and  v  »  T/2*  on  in  the 

usual  sense.  So  v  ■*  t/2»  on  in  the  sense  of  trace  of  lemma  8.1,  and 

v  e  K2 . 

Thus,  K2  is  a  closed,  convex,  nonempty  set  of  Hp{ft)»  and  a(v* ,v") 
is  a  continuous,  coercive,  real  bilinear  form  on  H^tft)  *  Hp(ft)  (Cryer  [I960, 
p.  54 9])  i  henoe  the  problem  (F2)  has  unique  solution  (Stampacchia  [1964]). 


Q  .E  «D  » 


Theorem  8.3.  If  k  <  kQ  then  the  problem  (F)  has  no  solution.  If  k  >  kg 
then  it  has  a  unique  solution. 

Proof :  If  k  <  kg  then  K2  is  empty.  By  proposition  7.2  K  is  also  empty, 

and  the  problem  (F)  has  no  solution . 

If  k  >  kQ,  then  take  v  as  in  (8.6).  We  have  known  that  v  e  Hp(ft) 

2 

and  v  »  T/2*  on  I*  .  But  from  (8.7)  we  see  that  |Vv|  <  kr  in  ft.  Hence 
v  e  K,  and  K  is  nonempty.  By  the  similar  argument  to  that  in  the  proof  of 
Proposition  8.2  we  obtain  that  the  problem  (F)  has  a  unique  solution. 

Q  .E  «D . 


Remark  8.1.  By  virtue  of  theorem  8.3  and  Proposition  7.2  we  obtain  that  the 
theorem  6.2  in  Cryer  [1980]  means  that  the  problem  (F)  and  (F,)  are  equivalent 
under  the  conditions  described  there. 
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Remark  8.2  .  The  conjecture  that  the  problems  (F)  and  (Fj)  are  equivalent  is 
not  right.  The  numerical  experiment  we  have  made  for  the  case  R(z)  =  1 
indicates  that  vQ  ?  v^  •  This  fact  may  be  shown  by  analytical  method  in  this 
rase  . 
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